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PREFACE 


This is the final "Final Report" for NASA Grant NAGW-690. This support began on 
December 1986 and ended on April 1995. Fortunately, limited funding from the Jet Propulsion 
Laboratory has allowed the continuation on a part time basis of our research efforts for NSCAT and 
ERS-1 until December 1995. 

This final report summarizes the research on model functions for ADEOS and ERS-1 and 
provides an analysis of the differences between three kinds of models. The success of the AMI on 
ERS-1 obtained at GSFC and NMC is highlighted. Predictions that were made prior to the receipt 
of an ESA Contractor Report are verified by that report. 

The problem of describing the wind stress is reviewed. 

A scatterometer model for high winds for the AMI on ERS-1 and ERS-2 is under 
development. The present status of this research is described. 

Finally, the accomplishments over the years of our past research are reviewed in terms of 
many published papers that were "firsts" toward the development of the methods used in theories 
that describe the winds, the waves and the development of the active microwave systems presently 
operational and planned for spacecraft. 

As this final report shows, there is much more that still needs to be done. How and when it 
will be done and who will do it remain to be decided. 
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THE PRESENT STATUS OF BACKSCATTER MODELS 


The Model Function Subcommittee 

A major effort during the past two years has been to contribute portions of the text of the 
document being prepared by the model function sub-committee which is charged with 
recommending the preliminary version of the "Model Function" {Long, 1995) to be used when 
NSCAT on ADEOS is placed into orbit. The portions of that report provided by us pointed out 
various problems that arose in the development of the model function for the SEASAT scatterometer 
(SASS). The document entitled "Current Progress in Ku-Band Model Functions" prepared by 
means of contributions from the model function subcommittee and dated 1 May 1 995 has been 
released along with a letter from the subcommittee chairman, David Long. It has been circulated 
to the full science team. Many of the contributions that we prepared have been included, but many 
others contributions were not. The report has many parts that need clarification and other parts that 
need corrections. Comments will be provide to Professor Long during the next few months. 

Many of these problems were documented by the research of Peter Woiceshyn in 
collaboration with Morton Wurtele. Since Mr. Woiceshyn was not a member of the model function 
committee, it was necessary to describe these results in considerable detail. For example, the results 
of Woiceshyn and Wurtele showed that there was a systematic concentration of wind directions near 
the upwind pointing direction for one of the antennas for some ranges of wind speed and incidence 
angles. For other ranges of incidence angles, the wind directions were concentrated at 45° between 
upwind and crosswind. Also for some wind speed and incidence angles there were systematic 
biasses in wind speed. These effects plus our own research indicated that picking a model to be used 
and an algorithm based on the maximum likelihood estimate plus assorted interpolation schemes 
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might be premature. The data from NSCAT should be analyzed very carefully to be sure that the 
best possible model is used before the data are released. 

The best possible model for NSCAT does not yet exist, and the present methods for 
developing one when data become available are not adequate. Unless great care is taken, the model 
will not be good for light winds and calms and winds over about 17 m/s. Our research, as described 
below, for the C-Band scatterometer on ERS-1 shows that it ought to be possible to obtain a model 
for winds as high as 27 m/s. 

Present C-Band Analytical Models 

The research done on the model function for the C-band scatterometer on ERS-1 illustrated 
both the success and the failure of the methods that have been used. In our opinion, none of the 
proposed model functions is optimum. The published papers on the models contradict each other 
and disagree substantially outside of the range of moderate wind speeds. Very high winds are not 
modeled correctly. The backscatter values for light winds are incorrect because of the way the data 
were processed before they were released. 

To illustrate the wide variation in the properties of the analytical functions that attempt to 
relate wind speed to backscatter for C-Band a simplified summary of some of them follows without 
naming the models. 

1) The first model is a power law model of the form 

log o°, = log A, + B, log V (1) 

or o° 1 =A, v® 1 (2) 
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for, say, upwind the graph of this model on a double log plot is a straight line for a fixed incidence 
angle. If the wind speed is zero, the backscatter is zero. 


2) The second model in simplified terms, for at least one version, can be represented by 

log o° 2 = log A 2 + B 2 Vv (3) 

or o° 2 = A 2 10 B ^ V (4) 

The backscatter equals A 2 if the wind speed is zero. 

3) The third model has no convenient analytical expression but its behavior relative to the 
other two can be described. 

o° 3 = o° 3 (V, V T ) (5) 

The backscatter equals zero if V = V T . 

All three of these models can be made to agree at two points when o° is plotted as a function 
of v. If then 

o°, (v,) = o° 2 (v,) = o° 3 (v,) (6) 

and o°j (v 2 ) = o° 2 (v 2 ) = o° 3 (v 2 ) for V, < V 2 (7) 

These three models then have the property that: 

For v<v,; o° 2 (v) > o°,(v) > o° 3 (v) (8) 

For v, < v < v 2 ; o° 2 ( v ) < o°,(v) < o° 3 (v) (9) 

For v 2 <v; o° 2 (v) > o°,(v) > o° 3 (v) (10) 

Each model can in principle be fitted over the limited range for v, to v 2 , say 5 to 15 m/s. 
However, if any one of the three should prove to be correct the other two will not yield good results 
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for low and high wind speeds. 

The difficulty with this approach for a model to relate backscatter to wind speed is that once 
the functional form is chosen, the process of fitting the function to the data cannot change that 
functional form and the results may fail outside of the range wherein most of the data lie. 

Figures 1 and 2 are schematic graphs for these three different models. Figure 1 is for linear 
scales. Since most data for winds fall in the 5 to 1 5 m/s range, it is very difficult to discriminate one 
model from another. For both higher and lower wind speeds, the three models diverge but 
nevertheless all three predict increasing backscatter with increasing wind speed. 

Figure 2 illustrates the difficulties of working in log space, especially for low winds. Model 
1 predicts a backscatter value of 0.01 for a wind near one millimeter per second. Model 2 predicts 
a backscatter value of about 0.012 for no wind. Model 3 predicts a backscatter value of 0.001 for 
a wind of 2 m/s. 

The Success of ERS-1 

The value of scatterometer data has been conclusively demonstrated by some unpublished 
results of Robert Atlas of GSFC and by the development of a method for the routine assimilation of 
AMI data from ERS- 1 into numerical analysis procedure at the National Meteorological Center. The 
data from ERS-1 are good enough to be useful on a day to day basis from a moderate range of wind 
speeds. 

During the NSCAT meeting in Hawaii, Robert Atlas showed the value of ERS-1 data in 
many ways. Low pressure centers were found that were not located by conventional data. Those 
that were located by conventional data were frequently not placed in the correct location. An initial 
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figure 2. DOUBLE LOG PLOT FOR THREE MODELS 





value update for a forecast can not yield a good forecast if it is incorrect. Atlas showed that the 
range of accuracy of a forecast could be extended by means of ERS-1 data. 

Gemmill, et al. (1994) developed a way to intercompare six different model functions for the 
AMI on ERS-1 . They were CMOD 4 from ECMWF, CMOD 51 from IFREMER, CMOD 5L from 
ESA (the European Space Agency) CMOD 6 from the University of Hamburg, CMOD 7 from 
NASA-JPL/OSU and ESA (CMOD 4) the ESA "fast delivery product". The model, CMOD 4, is 
a modification of the CMOD 4 used for the ESA fast delivery product. 

The six candidate model functions were compared in terms of objective scores by which the 
models can be ranked. This procedure does not necessarily mean that a better model cannot be 
found. In fact, it is already known that the one selected does not yield realistic values for high 
winds. There is the danger that in the rush to produce results those that are released may be 
misinterpreted. 

Tables 3b, 4b, 5b, 6, and 7b from Gemmill, et al. (1994) rank the statistics obtained by using 
these six models for a year of data and comparing their winds to data buoy measured winds. For 
example, for Table 3b, CMOD 51 had the lowest speed bias, ESA had the lowest RMS speed 
difference, CMOD 51 had the largest speed bias - and so on. 

The grand sum given above Table 3b shows that CMOD 4 and CMOD 51 are virtually in a 
tie and that the other contenders are far behind. 

On the basis of this study, the ERS-1 data are being routinely processed at the ocean products 
center of the National Meteorological Center as described by Peters, et al. (1994). One can 
confidently expect that an impact of these data will be demonstrated in due course. 
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TABLES FROM GEMMILL, ET AL . 1994. 


GRAND 

CMOD4 

CMOD5I 

CMOD5L 

CMOD6 

CMOD7 

ESA 

SUM 

51 

49 

109 

127 

125 

116 


Table 3b 

TRANSFER FUNCTION RANKINGS 


NDBC and TOGA buoys, High Seas, All Data 
Space Box: 0.5 degree, Time Window: +/- 3 hours 

Dates 93 09 09 - 94 09 09 



CMOD4 

CMOD5I 

CMOD5L 

CMOD6 

CMOD7 

ESA 

SPD BIAS 

2 

1 

6 

5 

4 

2 

SPD RMS 

2 

2 

6 

4 

5 

1 

SPD COR 

3 

1 

2 

6 

3 

3 

DIR RMS 

1 

2 

3 

4 

5 

6 

VECT CORR 

2 

1 

3 

3 

3 

6 

FOM 

1 

2 

3 

4 

5 

6 


SUM 11 9 23 26 25 24 


Table 4b 

TRANSFER FUNCTION RANKINGS 


NDBC Mid-latitude, High-Seas, Winter Data 
Space Box: 0.5 degree, Time Window: +/- 3 hours 

Date 93 11 01 - 94 04 31 

CMOD4 CMOD5I CMOD5L CMOD6 CMOD7 ESA 

SPD BIAS 
SPD RMS 
SPD CORR 
DIR RMS 
VECT CORR 
FOM 


• 1 . 3 5 

1 "I 5. 

2 ' 1 2 

3 1 1 

2 1 2 

1 2 3 


5 5 1 

4 5. . 1 

6 4 4 

4 5 6 

5 2 6 

3 5 6 


' SUM 10 9 


18 27 .26 24 
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Table 5b 

TRANSFER FUNCTION RANKINGS 


NDBC Mid-latitude, 

High-Seas 

, Summer Data 



Space Box: 

0.5 degree, Time ' 

Window: 

+ /- 3 

hours 


Date 93 09 

09 - 93 

10 31 and 

94 05 

01 - 94 

09 09 



CMOD4 

■ CMOD5I CMOD5L CMOD6 CMOD7 

ESA 

SPD BIAS 

2 

1 

5 

4 

4 

2 

SPD RMS 

1 

1 

5 

4 

6 

1 

SPD CORR 

2 

1 

2 

6 

3 

3 

DIR RMS 

1 

3 

1 

4 

5 

6 

VECT CORR 

1 

1 

4 

3 

4 

6 

FOM 

1 

1 

3 

3 

5 

6 


SUM 


8 

8 

20 

24 

27 

24 




Table 6b 






TRANSFER 

FUNCTION 

RANKINGS 


TOGA Tropical, 

High 

-Seas , 

Winter Data 



Space Box: 

0.5 

degree, Time Window: 

+ /- 

3 hours 


Date 93 11 

01 

- 94 

04 31 






CMOD4 

CMOD5I 

CMOD5L CMOD6 

CMOD7 

ESA 

SPD BIAS 


3 

1 

5 

6 

4 

2 

SPD RMS 


3- 

2 . 

5 

6 

4 

1 

SPD CORR 


3 

1 

3 

6 

2 

3 

DIR RMS 


1 

2 

3 

4 

5 

6 

VECT CORR 


2 

1 

4 

5 

3 

6 

FOM 


1 

2 

3 

5 

4 

6 

SUM 


13 

9 

23 

32 

22 

24 


Table 7b 

TRANSFER FUNCTION RANKINGS 


TOGA Tropical, 

High- Seas, 

Summer Data 



Space Box: 

0.5 

. degree. Time Window: 

+/ 

- 3 hours 


Date 93 09 

09 

- 93 

10 31 

and 94 05 

01 - 

94 09 09 

4 


CMOD4 

CMOD5I CMOD5L CMOD6 

CMOD7 

ESA 

SPD BIAS 


1 

4 

5 

5 

3 

1 

SPD RMS 


2 

3 

6 

5 

4 

1 

SPD CORR 


2 

2 

4 

6 

4 

1 

DIR RMS 


1 

2 

2 

4 

5 

6 

VECT CORR 


2 

1 

5 

3 

3 

6 

FOM 


1 

2 

3 

3 

6 

5 

SUM 


9 

14 

25 

26 

25 

SO 
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The ESA Contractor Report 

The report by Stoffelelen and Anderson (1995) reached us after most of the material in this 
final report had been written. The material in that report substantiates the statements made above 
and verified the need for a model for wind speeds over about 15 m/s. 

The notation, MLE, originally proposed by Pierson (1989a, 1989b) is used, as in Eqn. (1.4) 
of the report for K P = 0.05 approximately. Actually K P is a function of three parameters a, p and 
Y and can be very large for low winds. The parameters o, p and y have never been defined for 
ERS-1. 

Figure 3.13 in the ESA report reveals two of the difficulties with scatterometry. Both the 
ECMWF wind field analyses and the scatterometer winds do not identify calms separately. Both 
plots show no wind speeds above 21 m/s. Both features do not correspond to reality. 

Figures 3.16a and 3.16b compare the winds from CMOD4 to the ECMWF winds from a 
boundary layer model for the Braer Storm. Comparisons of boundary layer model winds, unless 
great care is taken, are found to be too low when compared to winds measured by a data buoy. Thus 
the wind speeds in both of these analyses are too low. Nevertheless, the CMOD4 winds are 
systematically lower than those for the boundary layer analysis. 

Similar results are shown for a storm that occurred on 14 March 1993 in figures 3.17a and 
3.17b. The ECMWF analysis shows winds of 50 knots (+/-) over a large area, whereas the CMOD4 
winds show 40 knots. Again both winds are probably too low. 

A graph of CMOD4 at upwind for various incidence angles would differ by modest amounts 
from the graph for curve 2 in Figure 1. To recover wind speeds of 25 to 30 m/s would require 
backscatter values that are much too high and that do not occur for these wind speed. 
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Our preliminary analysis of the Braer Storm attempts to fit power laws similarly to curve 1 
in Figure 1 for the wind speed range above about 1 5 to 20 m./s. 

WIND STRESS 

The nonlinear equations that describe turbulence in the planetary boundary layer are just as 
complicated as those that describe water wave properties and perhaps even more so. Attempts to 
parameterize the stress of the wind on waves and swell have been frustrated by the complexity of 
the problem. 

This was one of the objectives of our research during the past few years. No definitive 
solution has been obtained. The "Model Function Subcommittee" has recommended that the 
effective neutral wind at 10 meters be used because this value is not too sensitive to the assumed 
equation for the stress. The boundary layer models for computer based numerical weather prediction 
also do not seem to be very sensitive to the choice of a model for the stress. 

With the successful deployment of TOPEX/Poseidon, the need for a correct parameterization 
of the wind stress becomes important so that ocean currents can be modeled. 

Many of the efforts to parameterize wind stress in terms of the time averaged wind speed at 
10 m have determined the values of A and B in (U. 2 /U(10)) 2 = A + BU(10) = C DN with the result 
that the literature contains values for these two "constants" that differ be large amounts. The reason 
for these large differences may well be the lack of homosecidity, i.e., representative accurate values 
of the required data including a full range of wind speeds and wave conditions as defined by wave 
spectra. 
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Three important papers on this subject are those of Donelan (1982), Donelan (1990) and 
Blake (1991). Donelan ( 1 982) stratified wind stress values versus wind speed and wave conditions 
for short steep waves and longer less steep waves for the same wind speed ranges. The drag 
coefficient was higher for the same wind speed for the short steep waves than for the longer less 
steep waves. This showed that wave conditions had to be considered in a model for the wind stress. 

Donelan (1990) reviewed the subject and used dimensionless variables based on the 
Buckingham PI theorem to find a relations between waves, winds and the wind stress. 

Blake (1991) used three different sources of data and found that the data from Smith (1980, 
1988) could be used to relate the stress to both wind speed and wave height. A copy of this paper 
is appended to this report. 

An important result of Blake was to use drag coefficients proposed by each author on the data 
for that author and the data for the other two authors. The proposed drag coefficient did well for the 
data on which it was based, but not on the other two independent data sets. Thus the three data sets 
lack homosecidity and the data are mutually incompatible. 

Our opinion is that the results of Blake need to be pursued further and that they would 
provide a useful relationship for the wind stress as a function of wave height and the wind at 10 
meters. 

Advancing technology may make it possible to do a definitive experiment similar to the work 
of Smith. The directional wave spectrum, the atmospheric stratification, the actual wind stress, a 
time history of the variables and various other essential parameters will need to be measured at a 
number of different locations for a year or more as to obtain a homogeneous data base representative 
of the full range of variability of the quantities involved. 
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A SCATTEROMETER MODEL FOR HIGH WINDS 


Introduction 

The complete description of the variables that link o° with the vector wind at the surface 
of the ocean has yet to be determined. Each model function that now exists as described in simplified 
terms in a previous section tries to emphasize the suitability and the importance of particular 
variables that connect o° to the vector wind. Empirically, o° = f(V, x, 0, •••) where V is the wind 
speed, x the relative azimuth, 0 the incidence angle and the ... are all the other known and unknown 
variables that might be involved. 

Scientists must work towards establishing as complete as possible a list of these variables 
that when put in their right combination should produce the long awaited results that scatterometry 
has made boast of earlier. In a more practical sense, there is the need of scientists to understand 
better the basic functional form of the relationship between the electronically measured radar cross- 
section o° and the wind velocity in space and time. A simplified relationship between o° and V 
will introduce a "user-friendly" atmosphere among users of marine surface wind products. As part 
of a negative commentary, the following observation is noted. Operational forecasters are just 
beginning to make prognostications of modem weather phenomena with the use of scatterometry 
as a method of data enhancement and/or correction. One of the manifest functions of providing 
scatterometry winds is to help weather forecasters and ocean modelers in setting up the most 
accurate and complete data base to study geophysical systems. 

Operational remote sensing of real-time winds are just beginning. The data that are utilized 
do well over a moderate range of wind speeds, but numerous publications have shown that winds 
above about 15 m/s are not correctly described. 
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To obtain the surface wind at a height of 1 0 m above the ocean surface from scatterometry 
requires a set of logistics to be in place. These are a number of co-located o° values of the ocean 
surface from at least three viewing geometries and a model function that gives a correct mapping of 
the radar cross-section parameter to the field of motion. The wind vector is the by-product of this 
experiment. It is well-documented that if the magnitude portion of the vector can be obtained, then 
the directions can also be worked out by making use of varying radar viewing geometries and other 
available data, for example, isobaric and streamline analyses of existing weather maps. 

Given the state of the art in electronic engineering, scientists have placed substantial 
confidence in the backscattered o° values that are received and processed from polar orbiting 
satellites. Skill in data management and analysis help researchers to find colocated measurements 
of a 0 . The somewhat illusive model function continues to make the wind retrieval task both 
complex and difficult. Scientists are, however, encouraged by some partial results that are obtained 
from their work with many of the existing model functions. This approach using the method of 
synthesis has, unfortunately, placed scientists short of the mark. What is needed is an accuracy 
(within specified limits) that allows the wind output to be usable in operational and research 
situations. Thus there is evidence that existing model functions need refinement. As is 
documented, some of them work fairly- well in the mid-range (5-15 m/s) of wind speeds. The low 
and high wind regime ends are much more problematic. Here is the suggestion that there may be 
other parameters that are important in these parts of the wind spectrum, that are, until now, not 
included in the formulation. The role of sea-surface temperature at the low end, and spilling waves 
and breakers at the high end are parameters to be considered. Then there is the question of the 
limitation of the models themselves. 
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Our research goals over the past few years have been directed along rather practical lines. 
Some of the short-comings that have been cited earlier in this report are being addressed. Efforts 
were concentrated on evaluating and comparing model functions. (See Sylvester 1 984, Sylvester 
et al. 1989, 1990). Studies of this nature provide an understanding of the pros and cons of certain 
model functions. 

Most of the evaluation techniques are related to self-consistency. Models should be verified 
against conventionally obtained winds that define the range of speeds obtained over the ocean. For 
the past two years, efforts have been centered on providing a suitable network of conventionally 
measured surface winds to be used as a basis for verification, and in particular, winds in excess 
of 20 m/s. Products from ERS-1 should be verified against high winds so that the validity of model 
function can be attested. 

A further motivation for our approach was based on the following ideas. 

* It is known that nature produces high winds over the ocean but the ability to measure them 
is greatly reduced because many ships generally avoid areas of storms. 

* Absolute confidence cannot be placed on winds obtained by ships as pointed out by Pierson 

(1990). 

* Most of the data buoys that provide the best wind products are placed in selected near-shore 
areas and may not lie in the path of the storm which is the subject of the inquiry. 

The project IN SEARCH OF HIGH WINDS came into existence in 1993. As conceived, 
it was intended to solve the important problem of finding suitable winds against which model 
functions could be tested and also to verily the capability of the instrument on ERS-1 as to its wind 
measurement capability. As pointed out before, obtaining a vector wind field by conventional means 
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proved prohibitive. The next best method was to resort to the use of corrected synoptic weather maps 
from which a wind field could be established. 

Three NMC synoptic maps became available in the middle of 1993. The deepest and 
probably the most intense extratropical cyclone in the North Atlantic in this century has to be 
something worthy of study. Firstly, all additional data relating to the storm had to be obtained 
before reanalysis could be done and then the utilization of PBL models for obtaining the winds 
from the pressure field should follow. Finally, the problem of working with scatterometer data in 
obtaining a model function could be addressed. The description of the methodology of carrying 
out the analysis of a high wind case is detailed below. 

The research sought to provide a network for verifying winds obtained with remote sensors 
by using the correct pressure field and exploiting the known relationship between the pressure and 
the wind. 

Motivation 

1 . There are a minimum of errors in the reported pressure field in synoptic maps for the 
North Atlantic and elsewhere. 

2. Synoptic weather maps for the North Atlantic produced by different meteorological 
offices tend to have good agreement. Where differences exist, reanalysis by the 
researcher can quickly correct the situation, especially, if more data become 
available. 

3. A "best" pressure map can be produced by merging data sets and map sets. This 
final map holds the input for the calculation of the wind. 
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4 . 


The known relationship between pressure and wind and the effect of warm and cold 
air advection are then utilized. For example, the geostrophic and the gradient winds 
can be obtained from the pressure fields. 


Method and Material 

1 . Obtain synoptic analyses for NMC and ECMWF plus any other synoptic data made 
available after map time. 

2. Compare products, correct data where necessary, and merge analyses and so as to 
produce a "final" reanalyzed pressure field. 

3. Obtain sea surface temperature fields. In this case, climatology may be sufficient. 

4. Calculate wind from PBL models using pressure gradients and air-sea temperature 
differences as the main inputs. 

5. Blend the derived 1 0-meter winds with the existing buoy and corrected ship data 
to establish a vector wind field. 

6. Carry out the intercomparison between remotely sensed model winds from 
scatterometry and those obtained from analysis just described above. 


PBL Models 

The calculated winds are obtained by utilizing the models of Cardone (1969, 1978) and the 
University of Washington (UW) {Brown and Liu 1982). 
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a) A marine PBL model proposed by Cardone calculates the wind speed and in flow 
angle at anemometer height above the sea surface. In particular, the ratio of the 
anemometer level wind speed to the geostrophic wind speed and the angle between 
the surface wind direction and the surface isobars are calculated. This baroclonic, 
stability-dependent marine boundary layer model has been in operation since the 
1970's. It is utilized internationally to obtain winds over the ocean. Table 1 is an 
example of an output for latitude 50° N for the Cardone model. The table gives the 
wind speed (first row) and the inflow angle (second row) for an isobaric spacing of 
4 mb and air-sea temperature differences as exhibited therein. 

However, the following important features as illustrated by Pierson (1978) will be 
included in the model. These are: 

(i) the effect of a cold low in increasing cross-isobaric flow for the surface wind; 

(ii) the effect of a cold low in increasing the speed of the surface wind; 

(iii) the effect of warm air advection in decreasing cross-isobaric flow for the 
surface wind. 

b) The University of Washington (UW) PBL model uses sea level pressure and 
temperature fields as input. 

This is a similarity model that incorporates stratification, variable surface roughness and 
thermal wind within it. 

Products from both models will be intercompared. Greater confidence will be placed in areas 
where there is good agreement. Where the two models disagree markedly, dynamic criteria will be 
adopted to decide which model will be followed. The one chosen must establish confirmity with 
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TABLE 1. Wind Speed in Knots and Inflow Angle as a Function 
of Air Sea Temperature Difference and the Pressure Gradient. 
( Cardone Model ) . 

LAT = 50 ISOBAR SPACING 4MB DTX = 0 000000 DTY = 0.000000 


DtFINT 

» 

SPACE 

DEG. 

LAT 

-16 

-8 

-4 

-2 

0 

i 

2 

4 
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113 0 

B 

mm 

103.9 

1 
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mm 
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Wm 

B 


BB 
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17.4 

BB 
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30.5 

78.0 


72.8 
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meteorological continuity. 


The Braer Storm 

The Braer Storm is an example of a suitable data set. 

A well-developed cyclone or storm may persist for several days. The intensity and direction 
of the wind are diagrammed on isobaric maps. Strong winds are associated with large gradients, 
and from one synoptic map-time to another, the steadiness with which the wind blows can be noted. 
It is rather common to observe high winds with a fetch of more than 500 km for which speed and 
the "from-which" direction remain statistically constant. The confidence that these vectors are 
suitable candidates for testing the above stated hypothesis is greatly increased, especially when 
there are actual buoy or ship reports in the area that agree well with the isobaric analysis. 

The central pressures and isobaric gradients in extratropical cyclones are a direct measure 
of their intensity. Applying boundary layer physics to these known properties can result in realistic 
estimates of the associated vector wind. 

The Braer Storm of 9 to 1 1 January, 1993 is recorded as the deepest North Atlantic cyclone 
of this century. The storm developed as a complex wave on January 9 and moved north eastward 
between NW Scotland and Iceland on January 10. It satisfies the definition of a meteorological 
"bomb". The 1200 GMT NMC analysis shown is Figure 3 estimates a central pressure of about 914 
mb. As a further reference. Figure 4 shows the UK Met. Office surface analysis for the same 
synoptic time. Both maps exhibit good agreement and are mainly characterized by very strong 
gradients surrounding the low pressure of the storm. It is also important to point out that even at 
distances 10 to 12 degrees of latitude south from the center, there are winds measuring at least 50 
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figure 3. Final NMC Analysis: 1200Z Sunday 10, January 1993 
The Well-Developed Braer Storm straddling the North 
Atlantic showing a dual low system with central pressures 
of less than 964 mb. and 920mb. respectively. The inset is a 
portion of the ERS-1 swath for Orbit No. 7796, Ascending 
Node 146.7, Southbound. 
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FIGURE 4. 

Atlantic Surface Analysis at 1800 GMT, Jan 10, 1993 showing cyclone track 
and central pressure (mb) at 6-hour intervals. Isobars are drawn every 4 
mbar. Courtesy of Stephen Burt, Mortimer, Berkshire 
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knots. 


In general, data gathered from actual ships, data buoys and coastal stations showed winds 
in excess of 50 knots (25.7 m/s). Some reports were between 60 and 70 knots (30.84 and 35.98 m/s) 
with gusts to about 80 knots (41.12 m/s). There were also gusts in excess of 1 00 knots (5 1 .4 m/s) 


Office at particular stations. 

observations related to the Braer Storm in the North Atlantic. 


reported by the UK Met. 

The following are 

Pressure : 

Lowest Pressure (62° N 15° W) 
Buoy 44768 (60.8° N 1 9.4° W) 

Ship V2QC (63.2° N 17° W) 

Ocean Weather Ship Cumulus 
Area around 6 1 0 N 1 7° W 


912 - 915 mb 

Pressure 930.4 mb [18 GMT Jan. 10] 

Pressure fall of 47 mb in 1 2 hr. 

Pressure 934.5 mb [15 GMT Jan. 10] 

Pressure fall of 42 mb in 1 2 hr. 

Pressure fall in excess of 24 mb in 3 hr. 

Pressure fall of 60 mb in 1 2 hr. 


Winds : 

• Strongest winds were confined to within about 280 km of the storm center. 

• The geostrophic wind was estimated well in excess of 100 knots (5 1 .4 m/s). 

• Reported surface winds ranged from 50-70 knots (27.5 - 35.98 m/s). 

• A gust of 105 knots (53.97 m/s) was measured at Roma, North Scotland. 

The storm area experienced sustained high wind which are shown by the time series of the 

vector winds for three stations in Figure 5. In particular. Ocean Weather Ship Cumulus shows 

sustained wind in excess of 50 knots for the period 0300 to 1800 GMT on January 10. Figure 6 is 
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FIGURE 5 • Abbreviated Weather Reports for Two Ships and a Station at Thorshavn in the Faeroe 

Islands. A Flag Represents 50 Knots. A Full Feather is 10 Knots and a Half Feather 
is 5 Knots. The Numbers on the Upper Right are the Atmospheric Pressure. For Example, 
804 Represents 980.4 mb. The Upper Left is the Air Temperature. Thorshavn is at 
62° 01' N, 6 46' W. The Weather Codes, Where Shown, are From the WW Group in the Report 
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FIGURE 6. Examples of Measured Winds From Ships Associated with 
the Braer Storm. Wind Barbs are in kt. 
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a summary of measured winds for the synoptic times given by (times/date). 

As outlined earlier, products of the calculated winds from the pressure field can be mixed 
carefully with quality-controlled conventionally measured winds to produce a wind field against 
which satellite data can be compared. 

The high wind case has not been completely analyzed. On the one hand, the PBL models 
have been difficult to get. After securing particular versions of them, they needed to be 
reprogrammed and updated. On the other, the greatest disappointment came when, in searching 
through the ERS-1 tapes, it was found out that data over the Braer Storm was sketchy at best and at 
most times, non-existent. There were many large gaps in the o° data as the AMI was turned off at 
various points along the satellite track. However, the pass for Figure 3 has some data that can be 
used. 

Empirically Determined Models 

It is most important to have a better understanding of the basic functional form of the 
relationship between o° and the wind vector. One simple method is to plot a scatter of wind speed 
versus o° for a given incidence angle and wind direction. This has been tried often, but the results 
have not produced consistent results as shown by the many candidate models. There has been too 
much scatter in these plots and the regression equations have been inconsistent. 

In an area of steady wind, that is, an area in the satellite swath where wind direction and 
speed are invariant over all 19 cells in the 500 km. across-track region, a model function can be 
treated in ways to provide a simplified relationship between o° and the wind vector. The wind is 
allowed to change in the direction along track. 
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The Braer Storm provides a suitable candidates for this study. The assumptions 

a) that the wind does not change by very much in direction or speed across the width of the 
swath and 

b) that the wind only changes in the along track direction, should hold if meaningful results 
are to be obtained. 

As an example, Figure 7 describes a portion of the Braer Storm in which the marked area 
is covered by the footprint of AMI-ERS-1 orbit for Jan. 1 1 1993. As this is a northbound pass, the 
winds are blowing parallel to the aft antenna. In this configuration, value returned by the aft antenna 
are related to downwind and those by the fore antenna are at or near cross-wind. 

At each of the 19 cells in the cross track direction, the o° value and wind speeds are 
stratified according to incidence angle. This being accomplished the simple relationship associating 
o° with V can now be tested. 

AMI and o° Values 

The AMI is a 5 -cm active microwave instrument that measures the radar backscatter from 
the ocean surface. The capillary waves are the dominant transponders. 

AMI’s 3-antenna system is necessary to resolve the directional portion of the wind vector. 
A schematic of the antenna geometry is given in Figure 8 adopted from Gemmill et al, (1994), 
NOAA. This well-illustrated diagram highlights important information relating to sigma-0 and 
the antenna configuration. The figure depicts the antenna orientation for a northbound pass. A wind 
blowing parallel to, but opposite to the forward direction of an antenna is designated as upwind. 
When it blows just parallel to and in the same direction of advance of the antenna, it is labelled as 
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FIGURE 7 . NMC Final Surface Analysis 0000 GMT 11 January. 

Isobars Describe the Western Half of the Braer Storm. 
The inset is the Area Covered by the Footprints of AMI 
ERS-1 Orbit for 11 January 1993/ Asc. Node 330.7 Long. 
00:27 GMT. The Inset Also Defines the Study Area. 
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(Space Craft) 
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Only odd-numbered 
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block of data are 
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which is not drawn 
to scale. 
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centers of cells in 
consecutive blocks. 
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361 cells of scatterometer 
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blocks of wind data are 
separated by 25km. Grid 
spacing between cell center 
nodes = 25km. 


FIGURE 8 

Sketch of the ERS-I AMI 
Scatterometer anteiuia at sea 
surface, data geometery and 
organization'-- not to scale. 

9 = incidence angle at 
the surface cell site 
along the antenna trace. 
<J) = antenna- relative wind 
direction at the cell site. 
(Adopted From Gemmill, 
1994 ) . 


et al 
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downwind. Crosswind is near 90° to upwind or downwind. 

The methodology outlined in the section on Method and Material has been utilized to process 
the wind field for an area of the Braer Storm that coincided with the footprints of AMI (marked in 
Figure 7). The NMC synoptic map for January 1 1 at 0000Z shows strong winds north of 45°N. 

Actual ship reports in the vicinity of latitudes 45°N to 50°N are for 40 knots ( ~ 20 m/s). The 
gradients north of this area are tighter and therefore suggest the existence of winds in excess of 50 
knots (~ 25 m/s). The wind array on which calculations and conlusions are based in this report was 
arrived at after processing and reanalysis that incorporated the following information data bases and 
meteorological models and practices. 

• NMC final analysis 1 1 January, 0000Z. 

• Actual ship reports. 

• Information from the Braer Storm (UK Met. Office). 

• Cardone's methodology (Refer to Table 1). 

• ERS-1 AMI scatterometer wind and sigma-0 products (JPL). 

Figure 9 shows a histogram of JPL wind speeds for the study area. Attention is directed to 
the calculated mean of 13.8 m/s. When compared with data from the other sources, it is seen that 
the product from JPL underestimated the representative winds by 5 to 7 m/s. 

Winds of about 26 m/s were detected in the array. Wind magnitude at discrete points in the 
array served as the input values to the empirically determined model. 

Some interesting and informative plots emerged. As examples. Figures 10 and 11 are 
presented in this report as part of our preliminary results. The plots are for incidence angles 38° and 
42°. The top graphs are the regression lines predicted by Eqns. (1 1) and (12) that follow. From 
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Longitude (Deg.) No of obs 



AMI- ERS Footprints for Orbit Jan 11, 1993, Asc. Node 330.7 Long 00:27GMT 

(A Segment in Study Area) 



Latitude (Deg.) 
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these two lines an estimate for upwind can be made as given in Eqn. (13). 

The bottom Figures (10b, lib) highlight the relationship between sigma-0 (max) and sigma- 
0 (min). The slight scatter about the regresssion line is an indication that this empirical model may 
produce some useful results. This statement is made in light of similar plots made in earlier studies 
by Sylvester, et al. (1989, 1990). Those max-min plots show the awesome scatter of the data. The 
points scatter from an outer edge all the way down to the 45° line. In the present plots points are 
concentrated at the outer edges only. 

Although the results are somewhat premature, they are, however, encouraging. It appears 
from this study that the AMI scatterometer can detect winds that are in excess of 20 m/s. In this 
case winds are detected as high as 26 m/s. 

The derivation of an empirical model for high winds follows. 

The model proposed by Long (1985) with an additional cos 3x term is given by 
o°(v, x, 6) = U(1 + b, cos x + b 2 cos 2x + b 3 cos 3x )/(l + b, + b 2 + b 3 ) 
where 

log 10 U = 1.877 - 0.1466 cos 0 + 0.00105 0 2 + y • log 10 V 
and y = 0.09885 - 0.506 0 - 0.000406 0 2 
For upwind, this become 

y(0) 

o°(v,x,0) = U = F(0)V TV 

This model is a simple power law as described previously. The partial derivative of o° with 
respect to x is zero at % = 0 (upwind) and at x = 1 80° (downwind) so that the model is slowly 
varying with x f° r these two directions. 
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SIGMA-0 (MAX) dB SIGMA-0 (dB) 


SIGMA-0 (dB) VS. LOG(IO) V FOR THETA = 38 Deg. 
S(Down) = -24.591 + 11.918 * log(10) V 
S(Cross) = -39.822 + 21.544 * log(10) V 



FIGURE 10a. 

Regression Lines for 
Downwind and Crosswind 
for 0 = 38°. 




Downwind 

Crosswind 


SIGMA-0 (MAX) VS. SIGMA-0 (MIN) FOR THETA = 38 Deg. 
S(max) = -2.671+0.544 * S(min) 



FIGURE 10b. 

Max-Min. Plot for 
Colocated Sigma-0 
Values, for 0= 38°. 


theta = 38 deg. 
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SIGMA-0 (MAX) dB SIGMA-0 (dB) 


SIGMA-0 (dB) VS. LOG(IO) V FOR THETA = 42 Deg. 
S(Down) = -22.134 + 8.906 * log(10) V 
S(Cross) = -32.301 + 14.449 * log(10) V 



FIGURE Ha. 

Regression Lines for 
Downwind and Crosswind 


for 0 = 42°. 



Downwind 

Crosswind 


SIGMA-0 (MAX) VS. SIGMA-0 (MIN) FOR THETA = 42 Deg. 
S(max) = -3.379+0.527 * S(min) 



SIGMA-0 (MIN) dB 


FIGURE ii b . 
Max-Min. Plot for 
Colocated Sigma-0 
Values, q =r 42°. 



theta = 42 deg. 
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At downwind, Eqn. becomes 
v(0) 

o°(v, 1 80°, 0) = F(0) V rv 'R 
where 

R = (1 - b, + b 2 - b 3 ) /(I + b| + b 2 + b 3 ) 

For vertical polarization, the downwind upwind ratio is slightly less than one. 

Thus Braer Storm backscatter and wind speed data for the swath shown in Figure 7 for a line 
perpendicular to the track is very close to downwind for the aft antenna and close to crosswind for 
the fore antenna. 

A preliminary evaluation of the winds in this swath can be paired with backscatter values for 
each row of cells at constant 0 so as to obtain estimates of F(0) and y(0). The values for F(0) 
and y(0) given by Long (1985) yield results that are not consistent for the high winds in the Braer 
Storm so that it will be necessary to find new values for these two functions for each row of cells 
parallel to the subsatellite track. 

These figures show that the backscatter values in the data occur at the outer edge of the plots 
originally obtained from the SEAS AT data as in Sylvester, et a\. (1989, 1 990 ). This means that the 
downwind values and crosswind values can be estimated from the regression lines. 


Thus in base 1 0, 


log o c ° = G c + H c log V 

(11) 

log o D ° = G d + H d log V 

(12) 
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An estimate from upwind would be 


log 0u ° = Gu+ Hu log V (13) 

where G y would be slighly larger than G D , and would be slighly less than H D . 

The subscripts U, D and C refer to upwind, downwind and crosswind. 

As an extention of the results of Pierson given in the report by the model function 
subcommittee let, 

o B ° = log o° = (A 0 + A, cos x + A 2 cos 2 x + A 3 cos 3 x ) 

+ (B 0 +B, cos x + B 2 cos 2 x + B 3 cos 3 x ) log V (14) 

where o B ° is in bels and where A 3 = A,/3 and B 3 = B,/3 ensures that the minima of logo 0 will be 
at 90° and 270°. 

By requiring that Eqn. (14) reduces to Eqns. (11), (12) and (13), the result is 


Ao = (G,j + G d )/4 + G c /2 (15) 

A, = 3 (G 0 - G d )/8 (16) 

A 2 = (G u + G d )/4 - G c /2 (17) 

A 3 =(G u -G d )/8 (18) 

and B 0 = (Hu + H D )/4 + H c /2 ( 1 9) 

B, = 3 (Hu - H d )/8 (20) 

B 2 = (Hu + H d )/4 - H c /2 (21) 

B 3 = (Hu - H d )/8 (22) 

Eqns. (11), (12) and (13) can be converted to natural units as in 

o u "=« u V Pu (23) 
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o c °=a c V Pc (24) 

Od 0 =«dV Pd (25) 

Let 

o N ° =aV» (26) 

where 

a = o 0 + a 1 cos x + a 2 cos 2 x + a 3 cos 3 x (27) 

P = p 0 + p, cos x + P 2 cos 2 x + P 3 cos 3x (28) 

and 03=0/3 = p 3 = p,/3 (29) 

Then o 0 = (a y + a D )/4 + a c /2 (30) 

ai = 3 (ou - o d )/8 (31) 

“2 = («u + «d ) /4 - « c /2 (32) 

«*3 = (au - o d )/8 (33) 

Po = (Pu + PdV 4 + Pc/2 (j4) 

P, = 3 (Pu - p D )/8 (35) 

P 2 = (Pu + PdV 4 ' Pc/2 (36) 

P 3 = (Pu ' Pd)(8 (37) 


The graphs of these two forms for a model when plotted in either natural or logarithmic form 
will only be the same at % = 0°, X ~ 90° and 270°. For any value of the wind speed, for example 
for 


10 


log °B°_ „ o 

- On 


(38) 
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The graphs will be the same as V is varied at upwind, crosswind and downwind for any 0, but for 
other values of x they will differ. This property could yield improved values for the variation with 
wind direction by using either 

o B ° = log o° (combined) = K o B ° + ( 1 - K) log o N ° for 0 < K < 1 (39) 

on a logarithmic scale or as 

on = K ( 1 0° B ° + ( 1 + K )< for 0 < K < 1 (40) 

on a scale with natural units. 
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REVIEW OF OCEANOGRAPHIC RESEARCH AND NASA GRANTS 


The National Aeronautics Space Administration and the Jet Propulsion Laboratory have been 
the principal supporters of research at The City College since September 1973 along with a number 
of smaller grants and contracts from other sources. Prior to that, support was provided at New York 
University even before Skylab was put into orbit. 

Prior to that at NYU, Professors Gerhard Neumann and Willard J. Pierson Jr. were supported 
by the Office of Naval Research from about 1950 to 1965. 

It is pertinent in this final report to review the contributions to the subject of remote sensing 
that were made over the period of time from about 1949 to the present. Hundreds of scientists all 
over the world are using and improving upon the results first obtained by this research. 

However, research is always open ended. There are substantial areas where the present 
models and theories are inadequate. There are ways to improve on these models and theories. Much 
more still needs to be done. 

NASA has supported this research since 1965. The current series of grants from NASA 
began on December 1, 1986 for a total funding of $754,000.00 (NAGW-690). The present series 
of grants from the Jet Propulsion Laboratory began on August 20, 1986 for a total funding extended 
until December 31, 1995 of $673,999.00. The number of people supported at various times has 
varied. Many have moved on to do research at various NASA activities or as contractors to NASA 
and other sources of support. 

During the above time frame, this research has yielded a number of "firsts" all of which have 
had a profound effect on the present status of remote sensing. Some of them are listed 
chronologically below with a brief note on their impact to the subject. 
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Pierson, W. J. (1951): The interpretation of crossed orthogonals in wave refraction phenomena. 
Beach Erosion Board Tech. Memo. 21. The effects of wave refraction over the 
continental shelf made it difficult to hindcast waves for the Beach Erosion Board using 
the Sverdrup Munk method to calculate the significant wave height and period. 

Pierson, W. J. (195 1): On the propagation of waves from a model fetch at sea. Gravity Waves, NBS 
Circular 521 , 1 75-186. The first attempt to study the dispersion of waves. 

Pierson, W. J. (1952): A unified mathematical theory for the analysis, propagation, and refraction 
of storm-generated ocean surface waves. Parts I and II. College of Engineering, Res. Div., 
N.Y.U. Never published but in the gray literature. Linear superposition of Airy 
Waves. (See Kinsman 1965). 

Pierson, W. J. and W. Marks (1952): The power spectrum analysis of ocean- wave records. Trans. 
Amer. Geophys. Un., 33, 6, 834-844. This is the first paper ever published that showed 
how an ocean wave time history could be analyzed so as to resolve the total variance 
of the record into frequency bands. Spectra based on this result, as modified, are being 
routinely estimated from hundreds of wave records every day. 

St. Denis, M. and W. J. Pierson (1953): On the motions of ships in confused seas. Trans. Soc. Naval 
Architects and Marine Engineers, 61, 280-357. This paper revolutionized the study of 
ship motions in random waves. Papers that use this study appear routinely in the 
Naval Architecture Literature, but it is so well know that few of them refer to it as the 
source. The 1973 "Certificate of Appreciation" noted that the paper "Changed our 
thinking about the ship, the sea, and the ship upon the sea". Important papers on bow 
submergence and slamming, ship routing and data analysis of ship trials resulted from 
this paper. The theory is "linear". Efforts to do the nonlinear problem are presently 
under way. 

Pierson, W. J. (1955): Wind-generated gravity waves. Advances in Geophysics, 2, 93-178. 
Academic Press. New York. Defined waves as a stochastic process that still satisfied the 
linearized wave equations. 

Pierson, W. J., G. Neumann and R. W. James (1955): Practical methods for observing and 
forecasting ocean waves by means of wave spectra and statistics. H. O. Pub. 603, U.S. Navy 
Hydrographic Office, Washington, DC. The first attempt to forecast wave spectra. The 
book is still used at the United States Naval Academy. 

Pierson, W. J. (1960): (Editor and co-author) The directional spectrum of a wind-generated sea as 
determined from data obtained by the Stereo Wave Observation Project, by Cote et al., 
Meteor. Papers, 2, 6, New York University Press, N.Y. The first directional wave 
spectrum based on a linear model. See the books by Kinsman (1965) and Neumann 
and Pierson (1965). 
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Pierson, W. J. and J. F. Dalzell (1960): The apparent loss of coherency in vector Gaussian processes 
due to computational procedures with application to ship motions and random seas. Sci. 
Rept. under Contr. Nonr. 285(17) and 263(09), College of Engineering, Res., Div., N.Y.U. 

Introduced the concept of co and quadrature spectra to the study of ship motions. 

Pierson, W. J. (1964): The interpretation of wave spectrums in terms of the wind profile instead of 
the wind measured at a constant height. J. Geophys. Res., 69, 24, 5191-5203. The 

variation of wind with height needs to be considered in wave forecasting. 

Pierson, W. J. and L. Moskowitz (1964): A proposed spectral form for fully developed wind seas 
based on the similarity theory of S. A. Kitaigorodskii. J. Geophys. Res., 69, 24, 5181-5190. 

This proposed spectral form is still used to model fully developed seas. A paper in the 
1994 Transactions of the Society of Naval Architect and Marine Engineers found that 
this spectrum fitted very high waves data very well. 

Neumann, G. and W. J. Pierson (1966): Principles of Physical Oceanography . Prentice Hall, 
Englewood Cliffs, New Jersey. Provided the background on wave theory and wave 
forecasting up to the time of publication. 

Pierson, W. J., L. J. Tick and L. Baer (1966): Computer based procedures for preparing global wave 
forecasts and wind field analyses capable of using wind data obtained by a spacecraft. Sixth 
Naval Hydrodynamics Symps., Sept. 28, Oct. 4, 1966. Described the numerical spectral 
wave forecasting model used at the Fleet Numerical Oceanographic Center for many 
years. Some changes have been made. 

Greenwood, J. A., A. Nathan, W. J. Pierson, G. Neumann, F. C. Jackson and T. Pease (1969a): 
Radar altimetry from a spacecraft and its potential applications to geodesy. Remote Sensing 
of Environment, 1, 59-70. The first paper to propose radar altimetry. Compared the 
widely different results for the geoid over the oceans and showed how a radar altimeter 
could measure the geoid. 

Greenwood, J. A., A. Nathan et al., (1969b): Oceanographic applications of radar altimetry from 
spacecraft. Remote Sensing of Environment, 1, 71-80. Showed that a radar altimeter 
could measure the slope of the ocean across the Gulf Stream, tides, storm surges and 
tsunamis. All of the predictions of these two papers have been verified far beyond 
those that were made at that time with the exception of measuring tsunamis. 
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